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Abstract 

Wc introduce the notion of an approximation system as a generalization of Taylor ap- 
proximation, and we give some first examples. Next we develop the general theory, including 
error bounds and a sufficient criterion for convergence. More examples follow. Prerequisites 
are mostly elementary complex analysis. 

1 Introduction 

Analytic functions are characterized by the property that they can locally be approximated by 
a power series, viz. the Taylor series. When we are given an analytic function g on a domain U, 
and we denote its derivatives at a point xq G U by at := pW(io), then g is in fact the unique 
analytic solution for go in the following infinite system of differential equations: 

SiOo) = <H, 9i(x) = 9i+i(x) (x G U), 

where we assume all gi to be analytic on U. We could truncate this infinite set of equations 
after n steps, so that we get: 

gi(x ) = ai (i<n), g[(x) = g i+ i(x) (x G U, i < n). (1.1) 

In terms of integrals this can be restated as: 

rx 

g n (xo) = a n , gi(x) = a» + / gi+i(t) dt (x G U, i < n). 

To solve this system of equations, we can pick any analytic function for g n such that g n (xo) = a>n, 
and each choice leads after n integrations to a solution for g$. In particular, when we pick g n to 
be constant a n , the resulting solution for go will be the n-th order Taylor polynomial. 

In this paper we will introduce a generalization of the Taylor approximation, based on the 
following alteration of (jl.lj) : 

gi(x ) = at, g'i(x) = fi(g i+1 (x),x), 

where each fa is an analytic function in two variables. To get the generalized equivalent for 
the n-th order Taylor approximation, we truncate again these equations after n steps, so that 
rewriting them in terms of integrals gives: 

px 

g n (xo) ■= a n , gi(x):=ai+j fi(g i+1 (t),t) dt (xeU,i<n). (1.2) 

Jxo 
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By picking the constant function a n for g n , we get after n integrations a solution for go . Like 
was the case for Taylor approximations, the solution to the truncated sequence of equations 
(|1.2p serves as an approximation for an analytic function g which satisfies the complete set of 
differential equations As will be shown in Example 12.71 we get the Taylor approximation 

by simply putting fi(y,x) = y. 

The introduced approximation method is a generalization of Taylor approximation, in the 
sense that it allows for a possibly non-linear function /j at each single integration step. It can 
also be regarded as a generalization of Picard iteration, in the sense that if we take at = a and 
fi = f, i.e., if we take them identical for each i > 0, then essentially reduces to a single 
differential equation. The ra-th order approximation to it then coincides with the n-th Picard 
iteration, as will be shown in more detail in Example 12.81 

The first ideas about approximation systems and many of the basic properties and the exam- 
ples were already formulated by the first author in his BSc Thesis [3J. He introduces a related 
but different concept called expansion system in his MSc Thesis [5]. 

The contents of this paper are as follows. In Section [2] we define approximation systems as a 
framework for approximating analytic functions, and we illustrate their use with some examples. 
In Section [3] some propositions are proved about the effect of certain coordinate transformations 
on an approximation system. These will be useful in later sections. Next we prove in Section 
[JJ a general result on the error of the approximations. Section [5] starts with a discussion of 
the relationship between the coefficients of an approximation system and the derivatives of the 
approximated function at its basepoint. Next we use this relationship for the development of a 
criterion that will guarantee the convergence of an approximation system. Section [6] gives some 
more examples of approximation systems, including estimates of the maximal error. We end the 
paper with some concluding remarks. 

Notation 

B(a,R) denotes the open disk with radius R around a in C. By cj)° n is meant the n-fold 
composition of the function cf>. So 4>°° := Id and 0°( n+1 ) := <po(p on for n > 0. We write D\f and 
-D2/ for differentiation of a function / of two real variables with respect to its first and second 
argument, respectively. We write ||/||y denotes the sup norm of a function / on a domain Y . 

2 Definition and first examples 

Let U C C be a simply connected open set and fix xq € U (to be called a basepoint of U). For 
some given analytic function g on U we will define a new way of approximating g by means 
of a so-called approximation system. This system will involve a sequence {ai}[ =0 of complex 
constants, and a sequence {/«}[rg of analytic functions in two variables fa: Yi — > C, where 
Yi = Vi x U for some open set Vi C C such that a^i 6 Vi for i < r, and where fi(.,x) is 
non-constant for all x € U. In this setting, r can be a nonnegative integer, but we will also allow 
the possibility of r being infinite, in which case our sequences must be interpreted as {aj}^ 
and {/,}- . 
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Definition 2.1. The pair ({oj}[ =0 , {/i}£To) is called an approximation system (shortly AS) 
of order r in the basepoint xq for the function g if there exists a sequence {gj}[ =0 of analytic 
functions on U with go := g such that gi+i(U) C for all i < r and the following equations 
hold: 

0i(zo) = Oi (i<r + l), = fi(g i+1 (x),x) (xeU,i<r). (2.1) 

Then we will also say that the pair ({aj}[ =0 , {/i}£T ) is an AS for the sequence {<?i}[ =0 • 

Remark 2.2. Note that according to the definition of an AS, the pair ({aj}[ =0 , {/jj-^Q 1 ) de- 
scribes a set of differential equations which is satisfied by the sequence of functions {gi}[ =0 - 
Unlike most common systems of differential equations, there aren't any further relations im- 
posed between functions gi and g.j for different i and j. The whole set of equations, as described 
in (|2.ip . can also be formulated as 

9i(x) =a,+ ! fi(g i+1 (t),t) dt (2.2) 

for x G U and i < r. Now suppose that the pair ({aj}[ =0 , {/i}[=o) is given as in Definition 
12.11 but that functions gi and 5 are not yet given. Let for some n (0 < n < r + 1) an analytic 
function g n : U — > y„_i be given with g n (xo) = a n . Then unique analytic functions (z < n) 
exist on a sufficiently small simply connected open neighbourhood Uq C U of xo such that the 
pair ({oi}" =0 , {/ij^To 1 ) * s an AS ^ or 5 := 9® ^ or the sequence {<7i}™ =0 • The functions <?j can 
recursively be constructed by (|2.2p and we put 5 := go- 

Remark 2.3. If in Definition 12.11 only the second part {fiYiZo °f the approximation system is 
given, but not the aj, then we can still have a sequence {giYi=o satisfying the second equation 
in (|2.ip . Next the first equation in (|2.ip can be used as a definition for the ctj. With {aj}[ =0 
thus chosen the sequence {<?j}[ = o W1U satisfy the full Definition ^. II For this reason we will often 
omit explicit values of the dj, since they will be implied. 

Let us now apply the idea of Remark l2.2l bv choosing g n , which we will now call g^\ identically 
equal to a n . So let be given an AS as in Definition 12.11 fix n (0 < n < r + 1), and recursively 
define 

9 [ n\x) ■= a n , 

gf\x):=a t + f U (g^ (t) , t) dt (i < n) 

Jxo 

provided this makes sense because of the inclusions 

9 l £i(U)cVi (i<n) (2.4) 
being valid. Then put 

We conclude that the pair ({aj}" =0 , {/il^Tg 1 ) * s an approximation system of order n for g^ 
in the basepoint xq. In view of approximation theorems to be given later in this paper, the 
functions g^ n \ as n gets bigger, often turn out to be good approximations for a function g as in 
Definition 12.11 This motivates the following definition. 



(x G U), 



(2.3) 
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Definition 2.4. The AS ({a;}[ =0 , {/i}-~o) is called a proper approximation system of order n in 

the basepoint xq if condition (|2.4p holds for the functions gj n ' defined by (|2.3p . If it is a proper 
approximation system of every order n < r + 1, we call it just a proper approximation system 
(or shortly PAS). 

Because in many cases the basepoint of an AS either can be understood from context or 
has an abstract value xq, we will usually not mention it in the following. In case of possible 
confusion we will write g: (U,xq) — > C in order to emphasize that xq is the basepoint. 



Remark 2.5. In most cases, the domains can be chosen in such a way that condition (|2.4p is 
satisfied. In general, for given n, we can always shrink U to a smaller open simply connected 
neighborhood of xq such that (j2.4|) is satisfied. 

Example 2.6. Suppose we have an analytic function g: U — > C. Then every sequence of invert- 
ible analytic functions {fiYiZo in one variable (i.e. fi(y,x) only depends on y) with sufficiently 
large range can be used for an approximation system for g. The sequence {gi}\=a ls uniquely 
determined by g = g and g i+1 := f' 1 o g[. 

Example 2.7. Taylor approximation 

When in Example 12.61 we let fi(y,x) = y for all i > 0, this gives rise to an AS for the function 
g: U — > C, such that gi = (and therefore aj = g^\xo))- Because 

,* n » (x-x ) k ~ l 

ai+ ^ akjj— j—^ dt = l^ a k- 



(Jfe-i-1)! ^ K (k-i)\ 

x o k=i+l y 1 k=i y > 



and 

we get by induction that 



9 [ n ] (x) = a. 



(x - x ) k sr^ ( k ), x y 



fcl K ' k\ 

k=0 k=0 

which equals the n-th order Taylor Approximation. Obviously, the fi and as chosen in this 
Example give us a PAS of any order in the basepoint xq. 

Example 2.8. Ordinary differential equations 

Consider the following ordinary differential equation (ODE) on a domain U : 

g'(x) = f(g(x),x), (2.5) 

where / : V x U — > C is an analytic function. Suppose we have a solution g to this differential 
equation on the domain U such that g(xo) equals a given value a £ V at a basepoint xq € ?7 and 
such that <7(i7) C V. Now with /j := / and := a we obtain an AS for the function g by just 
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putting gi := g in (|2.ip . Now, for fixed n and for U a sufficiently small open simply connected 
neighbourhood of xq, the Picard iteration scheme (see for instance [H Ch. 1, §3]) 

gW(x) :=a, 

gt\x):=a+ f *)# (* < ") 

J XQ 

makes sense with gf^^U) C V (i < n). We recognize (|2.3p and (|2.4p specified for our Example. 
Thus our choice for the and aj gives us a PAS of any order n in the basepoint xo, but with 
U possibly dependent on n. In fact, the general theory of ordinary differential equations tells 
us that, for U a sufficiently small neigbourhood of xq, we get a PAS for any order and we have 
uniform convergence of g^ to g (see also Remark 14. 7p . 

Note that it follows from fj2.6[) by induction with respect to n — i that we can say the following 
about a special case of Definition 12.41 

If fi = f, <H = a for alii then gf ] = gf^ = ■ ■ ■ = = 9 [n ~ i] ■ (2-7) 

Example 2.9. Functional differential equations 

As a generalization of the ODE (j2.5[) consider the functional differential equation (FDE) 

g'(x) = /(jo^),4 (2.8) 

on a domain U, where / : V x U — > C is an analytic function and <j) : U — > U is an analytic 
endomorphism on U. Assume that g is a solution of this equation with a given value g(xo) = a 
at a basepoint xo € J7, and such that g{4>{U)) C V. Also let <fr on denote the ra-fold composition 
of the function (f). So 0°° : = Id and c/>°( n+1 ) := 4> o on for n > 0. Put 5i = 5 o c/> oi : U -»■ V. 
Choose V5 C y := ^ such that <?(c/>°( i+1 )([/)) C V*. Then we see that the functions 

fi(y,x) := (<t> oi )'(x)f(y,<l> oi (x)) (i > 0, (y,*) e V> x 17) 

yield an AS for the functions #o( = #)> <?iiS2> • • • • Indeed, (|2.1|) can be seen to hold as follows. If 

g[{x) = fi(gi + i(x),x) for certain i then 

g' i+1 (x) = <f>'(x)g'Mx)) = <p\x) fi(gi+i(<f>(x)),<f>(x)) 

= <jj(x) f{gi + 2(x)A o{i+l) (x)) 

= {^ +1) )\x) f{g t+2 (x),^ +1 Hx)) = f l+1 (g t+2 (x),x). 

Note that a, = g((p oi (x )). 

Suppose that moreover 4>(xq) = xq. Then dj = g(xo) = a for all i. Also suppose that the AS 
obtained above is proper. Then we will show that 

g W 0^ = g%ff-\ in particular gf ] = g^ o </>°\ (2.9) 

For the proof first observe that the second equation in (|2.3j) is equivalent to 

(gP)\x) = {^)'(x)f{gfl 1 {x)A° t (x)), gP(xo) = a. (2.10) 



(x e U). 



(2.6) 
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Now assume that o (j) = g^i^ \i m — j < n — i. (It is certainly true for m — j = 0.) Then 

(sf 1 o cf)\x) = 4>'(x) {gt ] )'{m) = <f>'(x) ^ oi )'{m) f{gt}Mx))^ o(i+1) (x)) 

= (^)\x)f( g ^(x\^^\x)) = (g} +1 [n + l])\x). 

Since also g| n ' o <fi and gj^} 1 ' have the same value a at xq, we conclude that (j2.9j) holds. Further- 
more, by (|2.10p we see that 

{g^)\x) = f(g^ 1 \ct>{x)lx). (2.11) 
Remark 2.10. Note that (|2.1ip together with (g^)(xo) = a is equivalent with 

{gW)(x) = a+ P f{g^m)),t)dt. 



■i'O 



Define an operator S acting on analytic functions h on U satisfying h(xo) = a and h(<p(U)) C V 
as follows: 



(Sh)(x):=a+ f(h(</>(t)),t)dt. (2.12) 

Then 

g W = S(g^) = S n (a). (2.13) 

Special cases of proper approximation systems coming from FDE's with xq being a fixpoint 
of 4> are given in Examples 16. 1| 16.21 16.41 and Remarks 16.31 16.51 There it is possible to generate 
the approximations g^ by (|2.13p using the operator S. Also Example 16.61 is based on a FDE, 
but as its 4> has no fixpoint, the results there cannot be obtained by iteration of S. 

It is interesting to compare our operator S given by ()2.12p and its iteration with the operator 
T occurring in Grimm [31 Proof of Theorem 1] . The FDE there is more complicated than (|2.8p 
(/ also depending on g' and <f> also depending on g), but our FDE can be obtained as a special 
case. For that case Grimm's operator T becomes 

(Th)(x) :=f(a + J^ x) h(t)dt,x). 

Then iterates T n (0) approximate g' rather than g and, for H{xq) = a, S and T are connected by 

(Sh)' = Tti, (S n h)' = T n ti, 

or in integral form: 

rx fx 

(Sh)(x)=a+ (Th')(t)dt, (S n h)(x) = a+ {T n h'){t)dt. 

Jx Jx 
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3 Coordinate transformations 

Proposition 3.1. Let ({aj}^ =0 , {fiYi=o ) ^ e an AS for the sequence of functions gi \ (U,xq) C 
(i < r + 1). T/ien a biholomorphic coordinate transformation <j): (U,xq) — > (U,xq) (mapping xq 
to xq) gives rise to a new AS ({ai} r i=0 , {fiYZo ) f or ^ e sequence of functions 'g i : (U ,xq) - >• C 
(i<r + l), where fi(y,x) := fc(y , <p(x)) <p' (x) andg^x) := gi o <f>(x) . 

Proof. Equation (|2.ip is satisfied because for i < r we have 

Si 0*0 = S'iO(») < /' / (» = fi{9i+l(<f>(%)), ^{x)W \x) = fi(g i+ i(x),<t>(x))(f)'(x) =J i {g i+1 {x),x). 

Moreover, g i has the same image as gi, which is contained in Vi, and f { has domain tyj x U , so 
that the constructed approximation system indeed satisfies the requirements. □ 

Remark 3.2. Let ({ai}[ =0 , {/i}i=o)> 5*' ^ ana - /i ^ e as ^ n Proposition 13.11 Assume that 
({ai} r i=Q , {fiYiZo ) is an ^-S for some sequence g i : (U,xq) — > C (i < r + 1), where the are 
not necessarily as in Proposition 13. II But if ~g n = g n o ^ for some n < r + 1 then by Remark 12.21 
and by the proof of Proposition 13 . 1 1 we still have gi o <j) = g t for all i < n. 

This Remark leads us immediately to the following Lemma, which will be used in the next 
section. 

Lemma 3.3. Let ({aj}[ =0 , {/i}£T ) be a PAS and ({cii}i =0 , {fiYZo) ^ e an AS for the sequences 
of functions gi: (U,xq) C and 'g i : (U,xq) — > C (i < r + 1), respectively, such that these are 
related through a biholomorphic coordinate transformation <f>: (U,xq) — > (U,xq) (i.e., ~g~i(x) = 
gi o <j)(x) and fi(y,x) := fi(y,4>{x)) 4>'{x) ). Then g^f 1 is well-defined and equal to gf 1 ^ o (p for all 
n < r and i < n and the second AS is also proper. 

It is also possible to transform the AS by a linear transformation of the first argument of the 
fi, as is formulated in the following Proposition. Its proof boils down to just a straightforward 
calculation like in the proof of Proposition 13.11 

Proposition 3.4. Let ({aj}[ =0 , {/i}£Zg ) be an AS for a sequence of functions gi'. (U,xq) C 
(i < r + 1). Then fi(y,x) := afi(a~ 1 (y — b),x) (a ^ 0) gives rise to an AS for the sequence of 
functions ~g~i : (U, xq) — > C (i < r + 1), where <?j(x) := agi(x) + b. If the former system is proper, 
then so is the second, and we have: ~g^ (x) = a g\ n ^ (x) + b. 

It is not possible to formulate a version of Proposition 13.41 involving a non-linear transforma- 
tion ifii of the first argument of fi, as is seen from the identity 

(V>i-l ogi)'(x) = ipl-iigiix)) fi(g i+1 (x),x) (i > 0). 

4 Error estimates of the approximation 

In the following propositions, we denote by T(v, w) the set of Ci-curves in our domain U con- 
necting u with v, and we define a metric djj on U by 

du(v,w) := inf £(7), 
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where £(7) denotes the length of the path 7. More generally, let ttipi) denote the length of the 
initial part of 7 € T(v,w), from the point v to the point t. We now prove the following simple 
lemma, which we will use in Theorem [ 



Lemma 4.1. Let 7 € T(v,w) be a C\-curve in U . Then, for integer n > 0, 



du(v,t) n \dt\ < 



n + 1 



Proof. We have / d v (v, t) n \dt\ < / Z t {lT\dt\ = / s n ds = u; . □ 
J 7 V 7 Jo n + 1 



Theorem 4.2. 

ir-l 



A. Lei ({oi}[ =0 ) {/«}[=o ) ^ e an AS for the sequence {gi}[ =0 9 := 9o, where gi has domain 
U and fi has domain Yi = Vi x U , and where we assume that V% is convex. Let this AS also be 
a PAS of order n < r + 1. Then 

|5(x)-^k^)l<Pl/o||y ---pl/™-l||^ 1 ||5n-<7n(xo)|| C / ^E^L (x £ U). (4.1) 

Here D\ denotes differentiation with respect to the first argument and ||/||y denotes the sup 
norm of a function f on a domain Y . 

B. Assume that moreover n < r. Then 



\g(x) - gW(x)\ < \\Dif \\Y ■ ■ ■ Pi/.-ilk., \\fn\\ 9n+li u)xu {x G U) - (42) 



du(x ,x) n + l 

X) ~ <T '{*)\ ^ WUOUYo ■ ■ ■ ||-^lJn-l||Y n _i \\Jn\\g n+1 (U)xU 

Proof. 

A. We will show by downward induction with respect to i, starting at i = n, that for i < n: 

\ 9i (x)-g [ ?\x)\ < \\D 1 f i \\ Yi ---\\D l f n ^\\ Yn _ 1 \\g n -g n (x )\\u ^2^1 ( x € u). (4.3) 

Then the case i = of (|4.3|> yields (|4.ip . Clearly, fj4.3f) holds for i = n because </J^ (x) = a n = 
9n{xo). Now suppose that for some i < n (|4.3p holds with i replaced by i + 1. Then for i we 
have the following string of (in)equalities: 



\9i(x) -gT\x)\ 



< inf / IMg^^-Mg^^Wdtl 



< mf / Pi/iikbi+iW-i/l+iC*)!^! 

< inf / H-Di/illy IPi/i+ilk+i • • • ||-Di/n-i||y„_i |bn -ffn(»o)||c/ ^ 

7Gr(x ,a;) J- 



(n-i- 1)! 
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The first equality is by substitution of (|2.2p and ()2.3j) . In the second inequality we used the 
convexity of Vi and in the third inequality the induction hypothesis. Furthermore, by Lemma f4.ll 



-yer(x ,x)J 1 (n — i — l)\ ■yer(x ,x) (n - i)\ (n — i)\ 

The last equality follows by continuity and the definition of djj. Combine this result with the 
earlier inequalities. Then we conclude that (|4.3|) holds for i unchanged, and thus the statement 
follows by induction. 

B. Now we show by downward induction with respect to i, starting at i = n, that for i < n: 

\gt(x) -gf ] {x)\ < \\Ihfi\W, ■■■ Pi/n-ilk^ \\fn\\ 9n+m xu ^Xl-ty. {x G U) - (44) 
Then the case i = of (|4.4p yields fj4.2|) . Clearly, ()4.4p holds for i = n because 



\9n(x) - g [ n ] (x) \ = \g n (x) -g n (x )\ 



fn(g n +i(t),t)dt 



■f() 



<\\fn\\g n+l (U)xud(x ,x). (4.5) 



The proof of the induction step is analogous to what we did in part A. □ 

Remark 4.3. Because the sets U, Vq, . . . ,V n are open, the sup norms at the right-hand side of 
(|4.1|) and (|4.2|) may possibly be infinite, in which case the theorem becomes a trivial statement. 
We will often use Theorem 14.21 in cases where the sets U, Vq , . . . , V n are bounded and where the 
analytic functions under consideration extend to analytic functions on open neighbourhoods of 
the closures of these sets. Then the sup norms are finite. 

Remark 4.4. Substitution of the inequality (|4.5p for \g n (x) — g n (xo) \ in (|4.ip would have given 
a weaker form of the inequality (|4.2p . with the denominator (n + 1)! replaced by n! . 

If the domain U is starlike with respect to xq then du(xo,x) = \x — xq\ for all x GU. Hence 
we have the following corollary to Theorem! 



Corollary 4.5. 

A. Assume the same conditions as in Theorem \4-%[ and further assume that U is starlike with 
respect to the basepoint xq € U. Then 

\g(x) - gW(x)\ < ||£>i/o||y ■ ■ ■ pi/n-llk-x ll«7n " ffn(^o)||f/ ^ G U). (4.6) 

B. If moreover n < r then 

\ 9 (x) - 9 [n] (x)\ < ||iPi/o||y • • • Pl/n-llk-i ||/n|| g „ +l(f 7)xC/ ("^"J), ( X G (4.7) 
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Example 4.6. Consider the PAS based on the Taylor series as described in Example 12.71 
There /j = id and Vi = C for all i. Furthermore, g n = and is the n-th order Taylor 
approximation of g. Then fj4.6|) and ()4.7p respectively give: 



(«) _ „(")r,.„MI„ |X ~ Xp| 
r 

n+l 



|^)-ff N (^)l<il5 (nJ -5 W (^o)||a 

n 



\g{x)-^{x)\<U^\\ u '* *°' . 

[n + 1)! 

These estimates of the remainder term coincide with familiar estimates in Taylor's Theorem. 
Thus it makes good sense to consider this g^ as an n-th order approximation of g which 
generalizes the n-th order Taylor approximation. 

Remark 4.7. For a PAS with r = oo it would be desirable to have that g^(x) — > g(x) pointwise 
or uniform on some neighhbourhood of xq as n — > oo. In the generality of Theorem 14.21 this 
cannot be concluded from (|4.ip or (|4.2p . However, things improve if we suppose moreover that 
all functions fi are the same function / on Yi = Y = V x U. Then (|4.2p yields: 

IffW -5 l J wl < ll/l|y (n+l)! ( 4 - 8 ) 

which implies uniform convergence on Y if / and D\f are bounded on Y . Of course, this is 
still under the assumption that we are dealing with a PAS. In particular, it is required that 
9i+i(U) C V for i < n. If moreover aj = o for all i then the requirement simplifies by (|2.7p to 
^[ m l (jj^ V for m < n. Then, if the AS is such that g% = g for all g and if we take for V an 
open disk around a and for U an open disk around xq such that / and D\f are bounded on 
V X U, then we can shrink U to a sufficiently small open disk around xq such that (|4.8p yields 
for all n that ^^(x) € V for x € £7. Of course, this is a classical argument which was used 
for the convergence proof of Picard iteration, see Example 12.81 The error estimate (|4.8p is also 
classical in the context of Picard iteration, see the last formula in [TJ Ch. 1, §3]. 

5 Convergence criterion 

Lemma 5.1. Let 

9i(x) = fi(9i+i(x),x). (5.1) 

Then, under assumption of sufficient differentiability of the functions involved, we have for n > 1 
that 

gl n+1 \x) = (D 1 f l )(g l+1 (x),x)g ( i l[(x) + (D% f t )(g l+1 (x), x) 

n m 

+ E Y.^ D T k h){9i + x{x),x)p m>k {g' i+1 {xl . . .,gl n + - m+1) (x)), (5.2) 



m=2 k=l 

where p m k is a polynomial with nonnegative coefficients. 
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Proof. This is by induction with respect to n, where the case n = 1 is: 

Si' 0*0 = (Difi)(gi+i{x),x) g' i+1 (x) + {D 2 fi){g i+ i{x),x). □ 
As an immediate corollary we obtain: 

Lemma 5.2. Assume that for given /, i/te derivative g[ depends on the function g^i ac- 
cording to (|5.ip . and assume sufficient differentiability of the functions involved. Then, in a 
given point xq, the values g'^xo), g'( (xq), . . . , g!- n+1 \xo) are completely determined by the values 
g i+1 (x ),g' i+1 (x ),..., gll\ (x ). 

Proposition 5.3. Let ({ai}£_ , {fi}lZo) be an AS. Then: 

1. Equations (|2.ip completely determine gf (xq) for i + j < r + 1 in terms of the AS. 

2. For fixed {/i}£T the derivatives g^(xo) (j = 0,...,r) are completely determined by 
{at}i=o- 

3. If r = co and if the AS is for the sequence {gi}^ , then the functions gi, and hence g, are 
uniquely determined by the AS. 

Proof. We only need to prove part 1. Let i + j < r + 1. By Lemma l5?2l we see that gi(xo) , g'^xo) , 
. ..,gy\x ) are completely determined by a i: g i+1 (x ), g' i+1 (x ), . . . , gY^[ 1 \x ). By iteration, 
9i(xo),g'i(xo), . . . ,9i(xo) are completely determined by Oi, a i+1 , . . . , a i+j -i, g i+j (x ), hence by 

a i, ■ ■ ■ > a i+j- n 

Corollary 5.4. Let us have a PAS ({oj}[ = 0) {/i}i=o) as * n Definition \2.4\ Then, for n < r + 1, 

gW®(x )=g®(x ) (j<n). (5.3) 
Also, for fixed {fi}™= > the derivatives g^ n ^\xo) (j = 0,...,n) are completely determined by 

Proof. Let n< r + 1. Then ({a;}™ =0 , {fi}^) is an AS both for {gi}f =0 and for {sj n] }™ =0 . Then 
the result follows from Proposition 15.31 □ 

Remark 5.5. If in the PAS (KK =0 , {/i}i=o) for the 

sequence of functions <?j we have that 
xo = 0, Of = 0, fi(y,x) is even in y and in x, and ^ is odd, then we can do better and we 
conclude that for n < \r + 1 

o M(;)( ) = aW(0) (j<2n). (5.4) 

Indeed, it then follows by induction from (|2.3p that all gj 7 ^ are odd. Next, by the proof of 
Proposition 15.31 (where the dependence on the a, now can be omitted since these are 0) we 
have for i + j < r + 1 and j odd that g •(()), 5^(0), . . . , (0) are completely determined 
by 

ffi+i(0)jSi-fi(0)> • ■ ■ ;Si+i ^(0)i an d hence by iteration by S^ + i^_ 1 ^(0), and hence completely 
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determined by the /j. A similar reasoning applies to the gf^ for j < 2n — 2i — 1. Thus (|5.4p 
follows. 

Similarly, if in the PAS ({aj}[ =0 , {/i}£Tg) for the sequence of functions gi we have that 
%o = 0, fi(y,x) is odd in x, and g\ is even, then we can again conclude that for n < \r + \ 

()5.4p holds (for j < 2n + 1). Indeed, it then follows by induction from (|2.3[) that all 
are even. Next, by the proof of Proposition 15.31 we have for i + j < r + 1 and j even 
that gi(0),gf'(0), . . . ,5 t - (0) are completely determined by aj,# i+ i(0), aj+^O), . . . ,5^ 2) (0), and 
hence by iteration by aj, Oi+i, . . . , a i + ij- A similar reasoning applies to the gf 1 for j <2n — 2i. 
Thus the claim follows. 

We will now generalize to the case of an AS the idea of majorizing a power series in absolute 
value by a power series with nonnegative coefficients. 

Definition 5.6. An AS ({aj}[ =0 , {fijlZo) is called positive if 

a,>0 (i<r + l), (Dfj^/iXtHu.so) >0 (i < r, k, I G Z> ). 

A positive AS ({ai}[ =0 , {/i}i=o) is said to dominate an AS ({oi}[ =0l {/i}[=o) if 

H<5i (i<r + l), |(^^/i)(a i+ i, x„) I < (^/^fei^o) (» < r, k, I G Z> ). 

Proposition 5.7. 1/ ({ai}[ =0 , {/i }i=o ) * s a positive AS for {gi}^ =0 then 

9i J \x )>0 (i < r + 1, j < r-i + 1). (5.5) 

If moreover r = oo i/ien g-^(xo) > for all i and j. 

Proof. Let i < r+1 and j < r — i+1. From Lemma f5.1l (in particular equation (|5.2|) ) and from the 
assumptions we see that g[ (x ), g\ 2) (x ),..., g\ j) ) x ) > if g' i+1 (x ) , g^fo) , . . . , g\ 3 ^ {x Q ) > 0. 
Then (|5,5p follows by iteration. □ 

Proposition 5.8. Let ({Si}[ =0 , {/i}[=o) be a positive AS for \jjiYi=Q which dominates an AS 

A. T/ien 

Is^^o)! <9i\x ) (i < r + 1, j < r-i + 1). (5.6) 

B. // moreover, for some n < r + 1, 



then 



l^(x )\<g^(x ) (j£Z >0 ) 



C. Assume moreover that ({aj}[ =0 , {/j}[ =0 ) is a PA5 o/ order n and that ({ai} r i=0 , {/i}^ =0 ) 
extends to a positive AS ({aj}^ , {/i}^ ) / or {ffil^o- 

< fl^o) (» < n, i G Z>„). (5.7) 
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Proof. 

A. Let i < r + 1 and j < r — i + 1 and have assumptions as in the statement of A. 
If \gi%(x )\ < gi%(x ) for k = 1, . . . ,j - 1 then it follows by (JOJ that 

j — 1 m 

771=1 fc = 

j — 1 m 

^ E £ (^^r^)(5m,xo)p m , fc (^ +1 (x ), . . . ,s£r%o)) 

m=l fc=0 

Then (|5.6p follows by iteration as in the proof of Proposition 15.71 

B. The proof is analogous to the proof of A. 

C. Let i < re. Then |gj*'(xo)| = l a «l < &i = gi(xo)- Furthermore, for j > we have |<?f (xo)| = 

where the last inequality follows from the last statement of Proposition 15.71 Now 
(|5.7p follows by using B. □ 

Theorem 5.9. Let ({Sj},^ , {fi}°Z ) be a positive AS for {gi}^Z which dominates a PAS 

({ai}™ , {/i} 4 =o) f° r {fi)i=o- Take R > such that x e U if \x - x \ < R. Then gW -> g as 
n — > oo, uniformly for \x — xo\ < R. 

Proof For \x — xq\ < R we have 

\g(x) - 5 N {x) | = | J2 9 {X0) f ^ (s - 

j=0 ■ 7 ' 



7. 

j=n+l j=n+l 



where we have used (|5.3|) . ()5.6[) and (|5.7p . Now use that the power series of g(x) around xq has 
positive coefficients and is convergent for x = xq + R. □ 

Proposition 5.10. Let ({aj}^ , {/j}^ ) ^ e a positive AS for {gi}^ suc ^ that for some R > 
we /lave [/ = B(x ,R) and B(a i+ i,g i+ i (x +r)) C V* /or a// r G (0,i?). T/ien ({a,}^ , {/i}~ ) 
is also a PAS for {g.;}^ . 

Proof. Fix re G Z>o and r G (0,-R). We will show by downward induction with respect to i, 
starting at i = re, that gf 1 (x) < gi (x) if < x — xo < r and \gf^ {x) — a« | < gi+\ (xo + r) — Of+i if 
|x — xo| < r. For i = re we have ffn n '(x) = a n and the assertions are clear. Suppose the assertions 
hold for i replaced by % + 1. Then 

of ] (x) = a i + f (t),t) eft. 
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Hence, for < x — xq < r we have 

9i {x) <ai+ I fi (g i+1 (t) ,t)dt = gi(x). 

J XQ 

Also, for | a; | < r we have 

px 

\gf\xo + x) - Oi| < gf\xQ + r) - Oj < / fi(gi+i(x + r), t) dt = ^(x + r) - a» □ 



Remark 5.11. It is possible to apply Proposition 15.101 to a specific category of FDE's from 
Example 12.91 With the notation and conventions from that example, assume that a > 0, 
4>(xq) = xq, and that eft and / have nonnegative derivatives (the function / in both its arguments) 
in respectively xq and (a, xo). Then, it follows by the composition and product formulas for 
power series, that fi(y,x) = (4>° l ) (x) f(y, cf) ol (x)) also has nonnegative derivatives in (o,xo). 
Hence the AS associated to our FDE is a positive AS. By Proposition 15.71 it then follows that 
the solution g of (|2.8p has nonnegative derivatives in xq. Now assume that U equals B(xq,R) 
for certain R > and that g and <f> are analytic on some open set containing the closure of 
B(xq, R). Also assume that for all positive r < R we have 

(/>(B(x ,r)) C B(x ,r). 

Letting Vi = B(a,g o 4> 1 (xq + R) — a), the requirements for Proposition 15.101 are then satisfied, 
and so the AS associated to our FDE is a PAS of all orders. Therefore, application of Corollary 
I43B yields 

Dn+l 

\g(x) - g [n \x)\ < y\\ v W 2 )'\\u ■ ■ ■ \W n )'\\u Pi/Ik • • • Pi/lly^ ||/lk 



(n + 1)! ' 

where Yi = B(a,g o 0*(xo + R) — a) x (ff l (U). Since 4> is an endomorphism on U we have 

(0 oi )'(x) = ( p , {r-\x))ct ) '{r- 2 {x))---cp'{x) < u% (x & u). 

This implies the more simple estimates 
\g(x)-gW(x)\ < \\4>'U nin+1) \\D 1 f\\ Yo ---\\D 1 f\\y n _ 1 ||/|| yB __ (5.8) 

1 n(n+l) R nJr ^ 
- w'\[B{x 0t R) (W D lf\\B(a,g(xo+B)-a)xB(x ,R)) n \\f\\B(a,g(x +R)-a))xB(x ,R) ^ + -g, • (5-9) 

From (|5.9p it is clear that g^(x) converges uniformly to g(x) on .6(0, R) as n — > co. 



14 



6 Examples 

In some of the examples below we will use the Chebyshev polynomials T p of the first kind and 
U p of the second kind (see §10.11, (2), (22), (23)]): 

T p (cose) := cos(p0), T p {x) = \p £ k ^~_ m (^) V ~ 2 \ 

U (cos 0) - sin ^ + ^) U (x) - T ( - 1)fc(P " fc)! (2x)^ 



We will also use the polynomial T+ which is defined as: 

b/2j ( 

~k\(p-2k) 



r+(x) := r*%(*x), t p +(x) = i P J2 Tu^S 



k=0 

We deduce that 

T^ 1- (sinhx) = Tp^" 1 sin(ix)) = i _p T p (sin(ix)) = i~ p T p (cos(^TT — ix)) = z~ p cos(ip7r — ipx), 
which equals cos(ipx) for p even and i _1 sin(ipx) for p odd. Hence 



(sinh 



x 



cosh(px) (p even), 
sinh(px) (p odd). 



9(0) 


= 1, 


g(x) 


= e x , 




9(0) 


= 0, 


g(x) 


= sinhx 


(p even) 


9(0) 


= 0, 


g(x) 


= sinx 


(p even) 


9(0) 


= 1, 


g(x) 


= coshx, 




9(0) 


= 1, 


g(x) 


= cos X. 





We start with some examples which are special cases of an AS determined by a FDE, as has 
been described in Example 12.91 The underlying FDE's with their solutions are 

g'( x ) = (g( x /p)) p , 

g'{x) = T p + (g(x/p)), g(0) = 0, g{x) = sinhx (p even), (6.1) 

g '(x) = (-l) p / 2 T p (g(x/p)), 3(0) =0, 9 (x) = sinx (p even), (6.2) 

g'(x) = sinh(x/p) U p -i (g(x/p)) , 
g '(x) = -sin(x/p)U p - 1 (g(x/p)), 

Note that in (|6.ip and (|6.2() with p odd the only analytic solution g of the FDE with g(0) = 
would be the function identically zero. 

Example 6.1. Let p € Z>i , R > 0, and consider the AS obtained by the choices f(y,x) := y p , 
<t>{x) := x/p, x := 0, g{x)~:= e x , U := B(0,R), V { := B(l, exp(R/p i+1 ) - 1) C 5(0, exp(i?/p m )) 
in Example 12.91 Note that eft leaves the basepoint invariant. Thus a, = a = g(0) = 1 for all i. 
Furthermore, 

gi(x) = exp(x/p 4 ) (x G U), fi(y,x) = p~ l y v ((y,x) e 7 4 x 17). 
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By induction we see that g^_ k is a polynomial of degree (p k — l)/(p — 1). Hence is a 
polynomial of degree (p n — 1). 

The AS satisfies the conditions of Remark 15.111 Hence, it is proper and by (|5.8p we get for 
\x\ < R the estimate 

\g(x) - g [n] {x)\ < ] IPi/lk • • • \\Dif\\ Yn ^ \\f\\ Yn 



(n + l)! 

< p-g«(n+i) TT p exp((p - 1)W +1 ) • exp(R/p n ) ■ - 
j=0 v 



nn+1 



e R 



p \n{n-l) ( n + 1 )| ' 

which converges to as n goes to infinity. For p = 2, the first four approximations are given by 
5 [0] (x) = 1, 5 W (_) = _ + _ , 5 M (_) = i + _ + ^ + g , 

roi , . X 2 X 3 7x 4 X 5 X 6 X 7 

q w (x) = 1 + x H 1 1 1 1 1 . 

y v 1 2 6 192 192 2304 64512 

Note that formula (|5.3p is indeed satisfied: the power series of g^ n \x) and g(x) = e x coincide up 
to the term with x n . 

Example 6.2. Let p G 2Z>i, and consider the AS obtained by the choices f(y,x) := T^~(y), 
(j){x) := x/p, x := 0, g(x)~.= sinhx, U := B(0,R), Vi := B(0,smh(R/p i+l )) in Example \2M 
Note that <j> leaves the basepoint invariant, so a% = for all i. Furthermore, 

gi (x) = sinh(x/p i ) (x € U), fi{y, x) = p _i T+(y) ((y, x)eViXU). 

The AS satisfies the conditions of Remark 15.111 Hence, it is proper and by (|5.8|) we get for 
la; I < R the estimate 



b(x)- 5 N (x) | < MWJfW \\D 1 f\\ Yo ...\\D 1 f\\ Yn _ 1 \\f\\ Yn 



R n+1 



(n + l)! 

un+l 

< p^-(n+D _____ ^(sinh^n+i)) JJ r+'(sinh(W +1 )) 



3=0 



^ n(re+1) 7-^7 cosh^/p^ 1 )) fT^^K 
F (n + l)! v IF 11 1 = 1 cosh( J R/pJ+ 1 ) 

# n+1 sinhtf cosh(iZ/p n+1 )) ^r 1 

i? n+1 e fi 
2p _n(n-l) ( n + l)| ' 



< 
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For p = 2, the first four approximations are given by: 

3 3 5 7 

fl r[°](x) = 0, g^(x) = x, g^(x)=x-\ , 5 [3l (x) = xH 1 h 



6 v ' 6 120 8064 

Here, even better than in Example 16. 11 g(x) and (which are both odd functions) have power 
series which agree up to the (zero) term of degree 2n. That this holds for all n > was explained 
in Remark 15.51 



Remark 6.3. By applying a coordinate transformation (ft: U —tU, (ft{x) = ix on the set U = 
B(0,R) of the previous example, we obtain by Proposition 13. II and Lemma 13.31 a new PAS with 
respect to the function sequence cjk(x) = g% o(ft(x) = smh{ix/p k ) = isin(x/p k ), where fk(y,x) = 
(ft'(x) fk{y, 4>{ x )) = W~ k Tp(y)- Next by applying Proposition 13.41 with a = —i and b = 0, we get 
a PAS with respect to (jk{x) = sm(x/p k ), where fk(y,x) = P~ k T£ (iy) = (— l) p / 2 p~ k T p (y) (p is 
still assumed to be even) . The error estimate of the induced approximations can again be given 
by Corollary 14.51 It turns out to be equal to the estimate given in Example 16.21 



Example 6.4. Consider the AS obtained by the choices f(y,x) := smh(x/p) U p —i(y) (p > 2), 



(ft(x) := x/p, x := 0, g(x) := cosh a;, U := B(0,R), Vj := B(l,cosh(R/p> +L )-l) in Example E3 
As before, (ft leaves the basepoint invariant, so aj = 1 for all j. The AS is a positive system 
(use [21 10.11(27), 10.9(3)]) and it satisfies the conditions of Remark 15.111 Hence by (|5.9p we 
have the estimate 

\g( X ) -g [n] (x)\ <^2™(™+ 1 ) (||Z?i/|| B( i, coshi? _ 1)xB(0ifl) ) n ||/||B(l,coBhfl-l)xB(0,fl) ( ra+1 )| 

(x 6 B(0,R)). 

Again by the positivity of the system we have 

Wf\\B(i,coshR-i)xB(p,R) = f(cosh(R/p),R) = smh(R/p) J7 p _i(cosh(i2/p)) = sinhi?, 
||-Di/||B(i,ooBhfl-i)xB(o J flD = (D x f)(cosh(R/p),R) = smh(R/p)U p _ 1 (cosh(R/p)). 
Also observe that 

sinhi U> ^ooaht) = £ ( ^M") = S -^± (pC oth(pi) - cothi) < (p _ 1} {t > ), 

p at \ sinhi / sinhi sinhi 

since an elementary analysis yields that p coth(pi) — coth t increases from to p — 1 as t runs 
from to oo. Hence 

sinh(R/p) U'(co S h(R/p)) < (p - 1) 



sinh(ii/p) 



Altogether, 



fl RsmhR / (p - l)RsmhR\ 



p|^(«+i)( n + 1)! V sinh( J R/p) 
which converges to as n goes to infinity. 
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Remark 6.5. Somewhat similarly as in Remark 16.31 we can apply a coordinate transformation 
(j): x i->- ix: U — > U to the set U = B(0,R) of the previous example. Then we obtain by 
Proposition 13.11 and Lemma 13.31 a new PAS with respect to the function sequence (jj(x) = 
g-j o(f){x) = cos(x/V), where fj(y,x) = 4>'{x) fj(y,(/>(x)) = -p~ j sax{x/p)U p -\{y). Again, the 
error estimate of the induced approximations is the same as the estimate given in Example 16.41 
For p = 2 the first four approximations g^ n \x) in powers of sin(2~ n_1 x) can be computed as 
follows. 

g [0] (x) = l, 5 W(x) = l-8sin 2 |, g® (x) = 1 - 32 sin 2 | + 160 sin 4 | - ^ sin 6 | , 

gW(x) = 1 - 128 sin 2 — + 2688 sin 4 — - 21504 sin 6 — + sin 8 — - 2326528 sin 10 — 

y K ' 16 16 16 3 16 15 16 

425984 12 x 1048576 14 x 

H sin sm — . 

3 16 21 16 

Similarly, we can expand g(x) = cosx in powers of sin(2 _n_1 x), by using that 

cos(2mt) = (-l) m cos(m7r - 2mt) = (-l) m r 2m (cos(±7r - t)) = (-l) m T 2m (sini), 

hence cosx = T 2 n+i (sin(2 _n_1 x)) (n > 0). 

Both g^ and g are even functions and their expansions in powers of sin(2~ n ~ 1 x) agree up to 
the (zero) term with power 2n + 1. This is true for general n by the second part of Remark 15. 51 
There the argument was given for expansions in powers of x. But then the result for expansions 
in powers of sin(2~ n_1 x) immediately follows. 

Example 6.6. Consider the AS obtained by the choices U = V = C, f(y,x) := y, <j)(x) := 
x + Q _1 logQ (a € (0,1)), xo := 0, g(x) := e ax in Example 12.91 It follows that fi(y,x) = y, 
gi(x) = a l e ax , cn = a % . The AS is proper. One easily verifies that 

n — i h h n h h 

a x , r „i . . v — v a x 



5 r(x) = a i ^— p, hence g^(x) 



k\ ' ^ k\ 

k=0 k=0 

Thus the n-th order approximation coincides with the n-th order Taylor approximation. Con- 
trary to the previous examples in this section we have here a special case of Example 12.91 where 
<p has no fixpoint and still everything can be worked out explicitly. 

Example 6.7. Let p G Z> 2 . Put A p>0 := 1 and X Pti := (p* - l)/(p i+1 - p 1 ) for i > 1. Let 
R £ (0,1). Take 

U := D(0, R), x := 0, g {x) := log (jZ^) , 9i{x) ■= (« ^ X )' 

fi(y,x) := X Pii y p . 

Then we have a positive AS with ai = 1 — i^o- Choose Vi := B(l, (1 — R)~ x p j — 1) for the domain 
of /i. Then application of Proposition 15. 101 yields after some computation that ({oj}^ > {fi}i^o) 
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is a PAS. So by Corollary 14. 5[ B we have for x £ U that 
|^)-3 N (^)l<ll/olk---H/n-llk- 1 ll/nlk ' 



p-1 , \p 



(n + 1)! 

< n pv< (( x - ^r vi ) K ' Vn ((! - «) -: 

p n ( 1 \ n +p~ 1 i^+i 



(n + 1)! 



< 



(p-l) n V 1 -^/ (n + 1)! 



where we used that \ v i < l/(p — 1) for i > 1. For p = 2 the first four approximations are given 
by 



„2 ^3 



(x) = 0, gW (x) = x, {x)=x+ X - + X - 

roi. . x 2 x 3 9x 4 3x 5 3x 6 9x 7 

ff M [ x )= x+ + + + + + 



2 3 64 64 256 7168 

Remark 6.8. In the previous example, clearly U cannot be enlarged so that it contains a disk 
D(0,R') where R' > 1, for the singularity at x = 1 cannot lie in the interior of U. However, 
when we allow U to have a non-circular shape, then numerical simulations indicate that the 
actual domain of convergence can become considerably bigger than the open disk with center 
and radius 1 (although x = 1 still remains a booundary point of course). The precise domain of 
convergence we have not yet established, but it seems to have a quite irregular shape. 

Remark 6.9. In all examples in this section, we have been able to calculate the approximations 

[nl 

provided by the algorithm explicitly. However, in many cases the approximation g\ is a poly- 
nomial whose degree is exponential in the number n — i of steps made. Because the computation 
time of each new step depends at least linearly on the input size of g\ , also the computational 
effort tends to increase exponentially as the order of the approximation increases. For practical 
purposes, this can become problematic when n gets large, even when making use of software 
like Mathematica. It is therefore worth noting that the algorithm suits itself well for numerical 
implementation. This can be achieved as follows: suppose we want a numerical approximation 
gv 1 ' of g^ (which in its turn is an approximation of g) along a path 7, going from xo to x € U. 
The first step is to partition 7 into a vector 7 : {0, 1, . . . , N} — s> Im 7, consisting of only N + 1 
points from the original path 7. Let := y(k + 1) — j(k) for k = 1, . . . , N . Now we construct 
gf ] : {0, . . . , N} -> C as follows: 

gt\k) = a n (k = 0,...,N), 

gf ] {Q) = ai (i = n-l,...,0), 

gf\k) = gf\k - 1) + h(gUi(k - l),7(fc " l))Ajfe_i (k = 1, . . . ,N, i = n - 1, . . . ,0). 

The vector g^ := g^ then represents a function with value g^(k) at the point j(k). This 
approximation of g^ will become better as the partition 7 of 7 becomes finer. Unlike exact 
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computation of g^ n \ computation of g^ carried out up to a certain numerical precision does 
not get significantly more complicated with each increment in n. Hence, provided that the 
coefficients aj are already given or computed, and that the functions /, to be applied do not 
substantially increase in computational complexity, the computational effort of this numerical 
implementation tends to be linear in both the depth n and the partition fineness N, rather 
than exponential. However, one should also take into account that exact computation gives 
in principle much more accurate results than the above algorithm. To get an increasingly 
accurate approximation of a target function g using the numerical algorithm, it requires beside 
an increment in n also a simultaneous increment in N, which makes a comparison between the 
approaches of course less straightforward. 

Another practical advantage of the numerical implementation comes up in cases when g^ 
cannot even be expressed in terms of known mathematical functions (cf. Appendix]). Then, 
numerical calculation is the only option. 

Concluding remarks 

A first remark is about the assumption of analyticity for gi and fi throughout the paper. This is 
done for convenience. Many results in this article could just as well be formulated for complex- 
valued functions with continuous derivatives of sufficiently high order. However, note that the 
coefficients only depend on the behavior of g at xq (just as with the the Taylor approxima- 
tion). Hence any changes to the function g outside a neighborhood of xq will not affect the 
approximations given by the algorithm. This makes the algorithm in particular suitable for 
analytic functions. 

Another assumption which could actually be relaxed, is that the domain U must be open 
in C. This assumption is made primarily for consistency in the conditions under which various 
theorems are proved. But not surprisingly, the algorithm can just as well be applied, for example, 
on a closed bounded interval in M for a C k function g on the interval with only left respectively 
right derivatives at the endpoints. If we would have worked with analytic functions on a not 
necessarily open set U then this should have been interpreted each time as that the function is 
analytic on some open neighbourhood of U. 

The examples we have discussed in this article mainly focus on approximation systems that 
arise from a single FDE, within the general framework of Example 12.91 For them we could 
do explicit symbolic computations and we could prove convergence of the approximations on a 
certain neighborhood of the basepoint xq. The algorithm certainly also works well for examples 
not coming from a FDE, as we know from numerical simulations. But the main problem there 
is that it is more difficult to prove specific error bounds, or even to prove convergence at all. In 
a somewhat different formulation, a couple of such examples are given in the appendix of [5]. 
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